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r ' Abstract. We examine the reduction process of a system of second-order ordinary differential 

-y ■ equations which is invariant under a Lie group action. With the aid of connection theory, we 
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curves of the original system can be reconstructed from the reduced dynamics. An illustrative 
example confirms the results. 
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O ■ 1 Introduction 

>'■ 

^ ■ This paper is concerned with second-order dynamical systems, or in other words systems 

H . of second-order ordinary differential equations, which admit a Lie group of symmetries; 

the question it deals with is how the symmetry group can be used to simplify the system 
(reduction), and how, knowing a solution of the simplified system one can find a solution 
of the original system (reconstruction) . 

This is not of course a new problem: reduction and reconstruction have been studied in 
a number of different contexts in geometry and dynamics. We mention in particular the 
following topics in which the object of interest is a second-order system: 

• the geodesies of a manifold with a Kaluza-Klein metric, and the Wong equations 

• Lagrange- Poincare equations and reduction by stages [Jf [T3]: 
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• non-Abelian Routh reduction [TO] : 

• Chaplygin systems [21 [5] . 

We mention these studies in order to emphasise the fact that we are engaged in this paper 
in something different from any of them. Each of the hsted studies deals with a special 
class of second-order system — for example, with systems of Euler-Lagrange equations, 
that is, equations derived by variational methods from a Lagrangian. We, by contrast, 
deal with systems of second-order equations pure and simple: we make no assumptions 
about how they are derived, and make no appeal to properties other than the property of 
being second-order and being invariant under a suitable symmetry group. (For the sake 
of clarity, we should perhaps remind the reader that by no means all second-order systems 
are of Euler-Lagrange type.) This has two important consequences. First, our results are 
more general than any of those obtained above: indeed, in some sense they must subsume 
the main features of the results of any of these particular studies. We shall show in fact 
how a system with symmetries may be reduced to a coupled pair of sets of equations, 
one of second order and one of first order. Second, our methods must be likewise more 
general: one cannot use variational methods, for example, if one is not dealing with a 
system of Euler-Lagrange equations. Of course one cannot say much without invoking 
some type of auxiliary machinery; but we use the minimum possible, just a connection. 

So far as we are aware, the only previous attempt at anything like this bare hands ap- 
proach is [3] , which deals with control systems of mechanical type, but with only a one- 
dimensional symmetry group; we deal with non-Abelian groups of arbitrary dimension. 

When we come to discuss the second-order case we shall represent a system of second- 
order ordinary differential equations by a special kind of vector field. There is therefore 
a lot to be said for starting off by considering the question for arbitrary vector fields, or 
in other words for first-order dynamical systems. 

Before proceeding, we must again make it clear exactly what we are attempting to do 
here. Perhaps the most studied example of an invariant first-order system is that of a 
Hamiltonian system with a symmetry group. The symmetry of such a system provides 
an equivariant moment map that is invariant under the Hamiltonian dynamics, and this 
feature plays an important role in the reduction process (this is discussed for example in 
[11] : a Lagrangian version is given in [8], together with an account of the reconstruction of 
the solutions in both the Lagrangian and Hamiltonian context). By contrast, we deal with 
vector fields pure and simple; since we have no Hamiltonian structure at our disposal, we 
cannot appeal to properties of momentum maps and so on. 

Suppose given a dynamical system, represented by a vector field X on a manifold M, 
which admits a Lie group G of symmetries. Suppose further that G acts freely and 
properly on M so that M is a principal bundle over a manifold B with group G; let 
vr*^ : M ^> B he the projection. Since X is invariant under the action of G there is a 
vector field X on B which is 7r*^-related to X; this is the reduced dynamical system. 

As a process of reduction, however, this is clearly incomplete in the sense that there is 
no way of reconstructing the original dynamical system from the reduced one; any two 



G-invariant dynamical systems on M which differ by a vr^^-vertical vector field (which is 
necessarily also G-invariant) have the same reduced dynamics. To see what is at stake, 
let us introduce coordinates (x*,a;") on M such that the orbits of G, or in other words 
the fibres of M -^ B, are given by x* = constant; the x* may therefore be regarded as 
coordinates on B. Let us denote by Ea a basis for the fundamental vector fields on M 
generated by the G-action; then 

^" = ''-^ 

for some non-singular matrix- valued function {K^g). Suppose further that we have at our 
disposal a distribution on M which is transverse to the fibres and G-invariant. Such a 
distribution will be spanned by vector fields Xi of the form 

d 

for certain functions A". We may then write 

X = Y'Xi + Z'^Ea = Y'-^ + (Z" - KyiEa- 

The necessary and sufficient conditions for X to be G-invariant are that Y^Xi and Z^-Ea 
are separately G-invariant. In particular, the F* are independent of the x", so that 
Y'^d/dx'^ may be regarded as a vector field on the base manifold B: this is the reduced 
dynamical system X, of course. The integral curves of X are solutions of the differential 
equations 

x' = Y\ x" = {Z'' - A^,x')Kl 

The equations of the first set define the integral curves of the reduced dynamical system. 
The remainder can in principle be used to reconstruct an integral curve of the original 
dynamical system from a known integral curve of the reduced one. 

This description of the process is somewhat disingenuous: a fibre-transverse G-invariant 
distribution on a principal G-bundle is of course just a connection, or more accurately a 
principal connection. This observation gives us the opportunity to describe the reduction 
in a coordinate-independent way, as is clearly desirable; when we do so, moreover, the 
reconstruction step acquires a more transparent geometrical interpretation than is ap- 
parent from the description above. Our basic contention is that the simplest additional 
machinery that is required to give a geometrically coherent account of the reduction and 
reconstruction of dynamical systems with symmetry is a connection; and we aim to show 
how these processes work in that context for second-order dynamical systems. 

A second-order dynamical system can be represented by a vector field F on the tangent 
bundle TM of a differentiable manifold M, of the form 



(9x" dv'^ 

where the f" are the fibre coordinates. Given a vector field X on M, let us denote 
by X*^ its complete, or tangent, lift to TM and by X^ its vertical lift. Then in terms 



of the structure described above in the first-order case, we may express a second-order 
differential equation field F on TP in the form 



V 



r = t;*Xf + v^E^ + D'XJ + D^E^ . 

It turns out that in order for F to be invariant under the action of G on TM induced 
from its action on M, each of the three components v°'E^, D°'E^ and f*Xf + D^XJ 
must be invariant. The last of these represents a second-order dynamical system, albeit 
in a generalized sense; the last two define the coupled first- and second-order equations 
which constitute the reduced system mentioned earlier. It is our aim to explain how this 
decomposition arises, with the help of connection theory; and to discuss the processes of 
reduction and reconstruction from this standpoint. 

In the following section we discuss the first-order case in greater detail. In Section 3 we deal 
with the connection theory required for the analysis of the reduction and reconstruction 
of second-order systems, and in Section 4 we carry out that analysis. Section 5 is devoted 
to consideration of an example. 

It will become apparent that a particular kind of Lie algebroid, the Atiyah algebroid of 
a principal bundle, plays an important role in the theory. In fact one can locate the case 
discussed here in a more general framework consisting of Lie algebroids and anchored 
vector bundles. Investigation of this aspect of the matter continues. 



2 First-order systems 

As before, we suppose that M is a manifold on which a Lie group G acts freely and properly 
to the right; we denote the action by ^"^ : G x M ^ M. Then tt^' : M -* M/G = B 
is a principal fibre bundle and tt*^ o ^^' = tx^' for all g E G. For ^ G g, the Lie algebra 
of G, we denote by C^m the fundamental vector field corresponding to ,^ G 0, that is, the 
infinitesimal generator of the 1-parameter group ip^puc) of transformations of M. 

The G-action on M can be extended to a G-action ip'^'^'' : G x TM —^ TM on the 
tangent manifold r : TM -^ M, given by {g,Vm) >— ^ Tmip^{vm), for m G M, Vm G T^M. 
This action equips TM with the structure of a principal fibre bundle over TM/G, with 
projection tt^*^. Then, of course, tt^*^ o ip'^'^' = tt'^^" and Tn'^'' o -i/;™^ = Ttt*^, where 
Ttt'^' : TM — > T{M/G); on the other hand, r o ■?/;J'*'^ = ip^^ We also have an action ip"^ 
on the space of sections of TM -^ M, that is, on A'(M), the space of vector fields on M, 
given by 

^^(X)(m) = ^J-(X(^-,(m))). 

A vector field X on M is G-invariant if for all g E G 

X{ipg{m)) = V^™(X(m)), or equivalently ipf{X) = X. 

If X is G-invariant then [^M) AT] = for all C ^ 0- If G is connected, as we shall generally 
assume to be the case, this is a sufficient as well as a necessary condition for invariance. 



The fundamental vector fields satisfy ipg{ih[) = (ad^-i ^)m, where ad is the adjoint action 
of G on g. 

For all m G M, tt™ induces an isomorphism TmM -^ {TM/G)^m(^^), the fibre of TM/G 
over 11^' {m) G M/G, and thus also an isomorphism TM -^ {ir^'yTM/G. As a conse- 
quence of this property there is a 1-1 correspondence between invariant vector fields on 
M and sections of the vector bundle r : TM/G -^ M/G (see e.g. [I2])- The vector bundle 
r has the structure of a Lie algebroid: the anchor map g : TM/G -^ T{M/G) is given 
by ["^l I— >■ T'K^'{v) (here and below [■] represents the equivalence class of the argument 
under G-equivalence, or in other words its G-orbit), which is independent of the choice of 
V G l^;] because of the property Tvr*^ o ip™ = Tii'^'; the bracket of two sections of TM/G 
is given by the bracket of the associated invariant vector fields. With this Lie algebroid 
structure r is called the Atiyah algebroid of the principal G-bundle vr^ [7]. 

The fibre-linear map Tvr*^ : TM — > T{M/G) is surjective on the fibres. The kernel of 
the induced map TM -* {■n'")*T{M/G) is isomorphic to the bundle M x g ^ M; the 
identification of M x g as a subbundle of TM is given by (?n, ^) ^^ ^m('^)- 

A connection on tt*^ is a right splitting 7 of the short exact sequence 

^ M X g ^ TM ^ {Tr^'YTiM/G) -^ (1) 

of vector bundles over M. The corresponding left splitting TM — ^ M x g will be denoted 
by uj. We will write w for its projection on g. The distinction between uj and vj can 
be made clear as follows. If we identify M x q with a subbundle of TM then uo may be 
thought of as a type (1, 1) tensor field on M; we have uj{^m) = C,m, while tu^C^m) = ^• 
Needless to say, both u and tu vanish on im(7). The map 7 may be thought of as the 
horizontal lift, the g- valued 1-form w as the connection form. 

If zu satisfies w{%Ij™v) = ad^-i zu{v) the connection is said to be principal. Equivalently, 
principal connections are right splittings 7 with the property 'jlipg'm, v) = ilj'^'''''j{m, v) for 
all V G Tj^M {jri){M / G) . The condition for the connection to be principal when expressed 
in terms of uj is simply that it is invariant under the G-action on TM, that is, that 

u O ^™ = ^™ o u. 

The manifold M x q comes equipped with the right action g t— > (^^^, adc,-i); we denote 
by g = {M X q)/G its quotient under this action. We remark that g is the (total space 
of) the vector bundle associated with the principal G-bundle vr*^ by the adjoint action 
of G on g; it is often called the adjoint bundle. When we take the quotient of the exact 
sequence ([T]) under the action of G we obtain the following short exact sequence of vector 
bundles over M/G: 

^ g ^ TM/G ^ T{M/G) -^ 0, (2) 

which is called the Atiyah sequence [7]. If 7 is a principal connection on vr*^ then 
TT'^'^''{'y{m,v)) is independent of the choice of ?7i G [m] = t{v) because of the invari- 
ance of 7; if we set 7(^1) = vr^*^ (7(771, 1;)) then 7 : T{M/G) -^ TM/G is well-defined 
and satisfies ^ o 7 = id, and is therefore a right splitting of the Atiyah sequence. This 
establishes a correspondence between principal connections on vr*^ and splittings of the 



Atiyah sequence, which is actually 1-1. If 7 is a right splitting of the Atiyah sequence, 
the corresponding left splitting will be denoted by cJ. 

If w is the connection form of a principal connection on ti^'' and X a G-invariant vector field 
on M then w{X) is a g- valued function on M which satisfies w{X) o ip^' = ad^-i w{X). 
So the map m \-^ {m, w{X){m)) E M x g is constant on the orbits of the G-action, and 
therefore defines a section of — > M/G. 

We next describe the reduction of a G-invariant vector field. As we pointed out earlier, 
a G-invariant vector field X can be identified with a section X of TM/G, given by 
X(7r*^(m)) = 7r^*^(X(m)). If 7 is a principal connection then X can in turn be decomposed 
into a vector field X = gX on M/G and a section ZJ{X) of g. The relation between X 
and X is X(7r*'^(m)) = T7r*^(X(m)); X is the reduced vector field of X. 

The decomposition may be described in a slightly different way. Given a connection 7, 
any X G A'(M) can be decomposed into its horizontal and vertical components with 
respect to 7; the horizontal component is determined by a section of {7t^")*T{M/G), and 
the vertical component can be identified with ui{X) and hence with the 0- valued function 
w{X). When X is G-invariant and 7 is a principal connection, the horizontal component 
is the horizontal lift of a vector field on M/G, namely X; and the section of that w{X) 
defines is just uJ{X). 

The following fact about the integral curves of an invariant vector field is well-known. 
Suppose that t \-^ c(t) is an integral curve of X, so that c = X o c. Then the curve 
t \-^ c{t) = 7r'^'{c{t)) in M/G is an integral curve of X, that is, 

t = Xoc. (3) 

Indeed, c = Tn'^'' o c = Tvr*^ o (X o c) = X o [n^" o c). In fact an integral curve c of X is 
completely determined by the underlying integral curve c of X and a curve t ^^ g{t) in 
G. To see this, note that there is a unique curve c^ in M, the horizontal lift of c through 
c(0), such that 

• Tf projects onto c 

• cT(0) = c(0) 

• the tangent c to cf is everywhere horizontal (so that 7f satisfies c = 7(0'^, c)). 

Then since c also projects onto c there is a curve t h^ g{t) G G, with g(0) = e (the identity 
element of G), such that c(t) = '^l''u)C^{t)- Now let 6 be the Maurer-Cartan form of G: 
then t 1-^ 9{g{t)) is a curve in g. By differentiating the equation c = ip^'cf we see that g 
must satisfy 

c=v^r((^(^))M°c^+^'). (4) 

Butc = Xoc = Xo iij^'c'^) = ^J"^(X o c^), whence 



{emMoc' + c' = Xoc\ 



The first term on the left-hand side (which when evaluated at t is the value at c^(t) of the 
fundamental vector field corresponding to 0{g{t)) G q) is vertical, the second horizontal, so 
this equation is simply the decomposition of X into its horizontal and vertical components, 
at any point of cT'. In particular, 

e{g) = w{Xol^)^ (5) 

the right-hand side being of course a curve in g. This is a differential equation for the 
curve g, and has a unique solution with specified initial value. (When G is a matrix 
group 6{g) = gg~^; and the equation 9{g) = ^, where t i— > ^(t) is a curve in q, can be 
written g = C,g, from which the assertion is obvious. See for example [15] for the general 
case.) Thus the curve g is uniquely determined by equation ([5]) and the initial condition 
g{0) = e. 

We can conclude the following. 

Proposition 1. Given a principal connection 7, one can reconstruct the integral curves 
of the G-invariant vector field X from those of the reduced vector field X. In order to 
carry out the reconstruction one needs to solve successively 

c = X{c) fore 

c = 7(0^,0) forcP 

e{g) = zuiXo'c') forg, 

to obtain finally the integral curve c = ip^'lp' of X . 

We now give some explicit expressions for the decomposition. We consider first the con- 
struction of a basis of vertical vector fields, that is, vector fields tangent to the orbits of 
the G-action. There are in fact two possible choices, at least locally, corresponding to 
what are sometimes called, as in [2], the 'moving basis' and the 'body-fixed basis'. The 
reference is to rigid body dynamics; the point is that the body-fixed basis is invariant. 

Let {Ea} be a basis for g, and C^^ the corresponding structure constants. The moving 
basis consists of the fundamental vector fields {Ea)M- These vector fields are not of course 
invariant: in fact for any fundamental vector field ^m, i^fi^Ai) = (ad^,-! C,)m, as we pointed 
out earlier. We will usually write Ea instead of {Ea)M for convenience. 

The definition of the body-fixed basis depends on a choice of local trivialization of vr" : 
M — ^ M/G. Let U C M/G be an open set over which M is locally trivial. The projection 
TT*^ is locally given by projection onto the first factor in f/ x G ^ f/, and the action by 
ipg'{x, h) = (x, kg). The maps 

Ea:U^{Mxg)/G\u by x^[(x,e),Ej 

will give a local basis for Sec(g) = Sec((M x g)/G) over U. These maps can be considered 
as sections of TM/G — > M/G by means of the identification 

Ea G Sec(s) ^^ Ea-.x^ 7T''{Ea{x, c)) G Sec (TM/G). 



Recall first that the injection M x g — * TM is given by (m, ^) i-^ C,m{'>tT')- Further, it is 
clear that the two elements ((x, e),^) and {{x,g),a.dg-i ^) of M x g belong to the same 
equivalence class in g. This is in perfect agreement with the above identification, since 

Now sections of TM/G can be lifted to invariant vector fields on M. For the above 
sections, the invariant vector fields are 

Ea ■■ {x,g) ^ (adg-i Ea)M{x,g) = ^™((Ea)M(x,e)). 

Then, indeed, ip^{Ea) = Ea- The corresponding basis for Sec(M x g) is given by the sec- 
tions {x,g) \-^ {{x,g),a.dg-i Ea). By contrast the fundamental vector fields Ea, identified 
as sections of M x g — >• M, are given by {x,g) i-^ {{x,g),Ea). The relation between the 
two sets of vector fields can be expressed as Ea{x, g) = A'^^{g)Eb{x, g) where {A^^{g)) is the 
matrix representing ad^-i with respect to the basis {Ea} of g. 

In fact the body-fixed basis {Ea} associated with a local trivialization {-k^')^^!] = U xG is 
obtained just by transferring to (7r")^^f/ the right-invariant vector fields on G associated 
with the basis {Ea} of g. The moving basis, on the other hand, corresponds to the 
left- invariant vector fields on G associated with the basis {Ea}. 

Let us take coordinates (x*, x") on M such that (x*) are coordinates on U, (x") co- 
ordinates on the fibre. Then there are 'action functions' (so-called in |2j) such that 
Ea = Ka{x'^)d / dx'^ . The relation [Ea, -E-fe] = CaiyEc leads to the property 

The invariance of the vector fields Ea can be expressed as 

We can use these differential equations as another way of constructing a body-fixed basis, 
as follows. We seek local vector fields {Ea}, given in terms of the moving basis {Ea} by 
Eb = ^b-Ea where (A^) is a locally defined non-singular matrix-valued function on M, 
which are G-invariant, which is to say that [Ea, Eb] = for all a and b. Thus, as above, 
the Al must satisfy 

Ea{At)+G:X = 0- (6) 

This is a system of linear partial differential equations for the unknowns A'^. The inte- 
grability conditions 

[Ea, Eb]{Ai) + CiEaiAl) - CtMK) = 

are identically satisfied by virtue of the Jacobi identity. The equations therefore have 
solutions locally on M, and a solution can be specified by choosing a local cross-section 
of the G action and specifying the value of (^1^) on it; the natural choice, which we make, 
is to take it to be the identity matrix. The A'^ will then be independent of the x*. 

8 



A simple calculation shows that 

where the A^ are the components of the matrix inverse to {A^). On the other hand, if we 
write [Eai Ei] = ^C^b^c then the coefficients (7^^ must be G-invariant, since everything 
else in the equation is. It follows that the value of C^^ along any fibre of (vr*^)"^?/ — > U 
is the same as its value on the section which determines the local trivialization, that 
is, where g = e; ii we take (A^) to be the identity there we obtain (7^^, = C^^,, that is, 
[Ea,Eb] = —C^f^Ec (as one would expect). 

We now consider the horizontal vector fields. We have at our disposal the local coordinate 
basis {d/dx'} of X{T{M/G)y, we put 

^. = 7(^1 GSec(TM/G). 

The sections {X,, Ea} form a basis of Sec(TM/G). They can be lifted to a basis {Xj, Ea} 
of X{M), consisting only of invariant sections. Then 

f d 

Xi{x,g) =7 \{x,g),— 

If we set 

^. = ^ - l\{x\ x^)Eb = ^ - i1{x\ x'^)Al{x'^)E,, 

then invariance of Xj amounts to 

For future use we calculate [Ea,Xi\ here also. We have 

[K,X,] = -X,{Al)A%. 

Now dA\/dx'^ = 0, and so 

X,{A:) = -^fAlEMl) = itCtjA^Ai = ^fC^Al, (7) 

so that 

[4,X,]=7.'C7:,4, (8) 

assuming as we may that [Ea, Eb] = ~CabEc- 
A vector field X on M can be written as X = Y^Xj + Y^Eb- If X is invariant then 

^^^ - n d ^^^ -{] 
dx'^ dx'^ 

If alternatively we set X = Y^Xj + Z^Ec, where Z'^ = A^Y^, then the second invariance 
condition becomes 

E4Z'^) + Cl^Z^ = 0. (9) 

An invariant vector field X projects onto the section X : (a;*) 1-^ Y^{x'^)Xj + Y"'{x^)Ea of 
TM/G and the vector field X : (x*) t-^ Y^{x^)d/dx^ on M/G. Finally, we have 

9 



• UJ{X) = VEa = Z'^Ea G TM- 

• uj{X) = Y-Ea G Sec(0); 

• zu{X) = Z-EaeC^iM,g). 

In fact, as we pointed out before, when X is invariant w{X) defines a section of g — > M/G. 
We now wish to explain how one can recognise a section of g in terms of coordinates. 
Recall that a section of g can be thought of as a function M — > g which is constant on the 
equivalence classes of the equivalence relation defining the associated bundle structure; 



that is, a g-valued function s on M such that s o ip^' 



adg-i s. Assuming as always that 



G is connected, we may equivalently write this condition as ^m(s) + [^, s] = for any 
^ G g, where the bracket is the Lie algebra bracket of g. We may express s as s = s^-Ea 
with respect to a basis {Ea} of g; in terms of the components s" of s the condition for s 
to define a section is 

E,{s'') + C^.s'' = 0. (10) 

This makes clear the significance of equation (jH]). 



3 Second-order diagrams and connections 

In this section we discuss the connection theory relevant to second-order dynamical sys- 
tems. Before we do so, however, it will be convenient to make some remarks about 
splittings of short exact sequences in general; these remarks will be useful later. 

If — i> ker / -^ A — > 5 — i> and -^ ker g ^ B ^ C ^ are two short exact sequences 
of vector bundles over the same manifold, then the sequence 

-^ keiig of)^A^-^C^O 

is also exact. Moreover the restriction of / to ker ((7 o /) gives rise to a fourth short exact 
sequence 

-^ ker / -^ ker((7 o /) ^ ker ^f ^ 0. 

In summary, we can draw the following commutative diagram: 



ker/ 



kcr(g o /) 



ker ( 



ker/ 



A 



B 



C 



C 
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The following facts are immediate. Suppose given splittings 'yi : B ^ A and 72 : C ^ A. 
If we set 73 = / o 72 : C -^ -B then 73 is also a splitting. If, in addition, 72(C) C 71 (-B) 
then 72 = 7i o 73. Furthermore, 71 restricts to a splitting ker^f -^ ker{g o /). Therefore, 
given 7i and 72, each element of A can be uniquely decomposed into three parts, one in 
C, one in ker/ and one in kerg. We will use this observation for our decomposition of 
second-order differential equation fields. 

We can now turn to the principal matter in hand. The actions of G on M and TM 
induce also a G-action on TTM: ^^™ : G x TTM -^ TTM, {g,X^) ^ T„^™(X^). As 
before, this means that there exists a principal fibre bundle structure ir'^'^'^" -. TTM -^ 
TTM/G with the properties tt^^"" o ^J^^'^ = n^^^' and Ttt^*^ o ^J^^'^ = T7^^*^ Again, 
the fibres of TTM and TTM/G are isomorphic, so TTM ~ {tt^^'YTTM/G. Therefore, 
the maps [Tvr^^^] : TTM/G -^ T{TM/G), |X] ^ Tn'^^'{X) and iTTvr*^] : TTM/G -^ 
TT{M/G), |X] \-^ TTti^\X) are well-defined and lead to the following commutative 
diagram: 



T{TM/G) 



TTM 




Tg 



TT{M/G) 



The above diagram contains bundles over TM^ TM/G and T{M/G). All of them are 
Lie algebroids. For example, TTM/G is the Atiyah algebroid of the manifold TM. In 
addition, all the maps in the diagram are Lie algebroid morphisms. 

First, we will consider the outer triangle consisting of the spaces TTM, T{TM/G) and 
TT[M/G). Since vr^*^ is a principal fibre bundle, the kernel oiTir'^^" can be identified with 
TM X 0, by means of the identification {v, ^) \-^ ^tm{v) = C]^(w). The two bases {Xj, Ea} 
and {Xi,Ea} of '^{M), based on the moving and the body-fixed basis, can be used to 
construct the basis {X^,E^,X^, E^} of X(TM), which we call the standard basis, and 
also the basis {X^,E^,X^,E'^}, which we call the mixed basis. Since C,tm = ^mj i^ is 
clear that the vector fields E^ span the vertical subbundle of the projection tt^*^. The 
advantage of the mixed basis over the standard basis is that the vector fields Xf , X^ and 
E^ are all invariant: 



[E'a, xn = [K, x,r = 0, [E-, X-] = [K, x,]^ = 

K, ^.1 = [Ea, E,Y = 0, [E^^, E-] = [E^, E,] 



- GabEc ■ 

So the vector fields Xf ,X/ and E^ can be projected to sections of TTM/G (by means 
and also to vector fields on TM/G (by means of Tir'^^"); the latter, denoted 



of TT^^*^ 
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by Xf , X^ and -E^, form a basis of X{TM/G). The following remark may be of some 
interest. Observe that TM/G — > M/G is a vector bundle; one can therefore define 
a vertical lift operation taking sections of TM/G -^ M/G to vertical vector fields on 
TM/G. The vector fields Xf and E^ on TM/G are in fact the vertical lifts of the 
sections Xi, Ea € Sec{TM/G); so we could write X^ = X^ and E^ = E^. 

The vector fields {E^,E^} span the vertical subbundle of the projection Ttt*^; so the 
kernel of TTtt'^' is isomorphic to TM x Tq ~ TM x g x g, the isomorphism being 
X"E^ + Z^'E^ ^ (X"E„, Z"E„). The vector fields {E^} span the kernel of the projection 
g; so ker Tg is isomorphic to TM x g, by Z"'E^ \-^ Z^-Ea- 

In this way we arrive at the following diagram of short exact sequences (taking into 



account the fact that {TTT^'yTT{M/G) 



TT^ 



'g*TT{M/G)y. 



TM X 



TM X0 



— TM X QX g 



TTM 



Tn^ 



TTtt" 



[TiT^YTTiM/G) 
{TT^'"')*g*TT{M/G) 



TM X 



{tt^''')*T{TM/G) 



{Tr^")*Tg 



(Tr^'"')*g*TT{M/G) 



There is a similar diagram for the spaces over M/G: 



T*B 



r^9 



T*(0X0) 



TTM/G 



ITtt- 



ITTtt" 



T{TM/G) 



Tq 



g*TT{M/G) 



Q*TT{M/G) 



Recall that r is the projection TM/G -^ M/G. The identification {TM x g)/^ ~ t*q is 
given explicitly by [vm,i] ^ (bml, H,^l)- 

We will use the basis {^, -£-[} for X{T{M/G)). The basic sections {^, ^} can also be 
used for bases of vector fields along the projection in the case of certain pull-back bundles, 
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tha^is_botl^for Sec((T7r")*TT(M/G')) and Sec{g*TT{M/G)). Finally, we will also use 
{Xf,Jq,E^} as a basis for Sec((7r™)*T(TM/G')). 

The two square commutative diagrams above will play the same role as the short exact 
sequences ([1]) and ([2]) in the first-order case. We will show that a principal connection 
on M induces splittings for all the short exact sequences in the squares. The idea is that 
connections of the M square which are G-invariant in the appropriate sense automatically 
give rise to connections for the M/G square and vice versa. We shall explicitly construct 
connections in both the middle horizontal and vertical sequences in the M square diagram, 
and use the general results at the beginning of this section to complete the task. 

The first induced connection lives on the middle horizontal line of the M square diagram; 
it is the so-called vertical lift of the principal connection on M. 

Proposition 2. Suppose given a connection on a principal G-bundle ti^' : M — > M/G, 
specified by its connection form w. The pull-hack t^vj of vj to TM (where t : TM -^ M 
is the tangent bundle projection) is the connection form of a principal connection on the 
principal G-bundle TM -^ TM/G. 

Proof. Clearly, t*w is a g-valued 1-form on TM. The action ip'^'^' of G on TM is r-related 
to the action ip'^'' on M. Moreover, the fundamental vector fields corresponding to the 
two actions are related by ^tm = (^m)° for any ^ G g; and in particular Tt{C,tm) = C,m- 
Thus 

r*ro(^rM) = ^(Tt{^tm)) = ^(6/) = C, 

while 

^™V*ro = T*^^'*w = adg-i T*zu, 

as required. D 

The connection defined by t*w is called the vertical lift of the original connection; its 
right and left splittings are denoted by 7^ and u^ (so that the connection form w^ is just 
given by to^ = t*w). The right splitting 7^ at the level of the M square can be given as 
follows. Let 7r^^^(f ) = -u; then 

7^ : {t^^'')*T{TM/G) -^ TTM, {v, X^^) ^ W, 

where W is determined by the condition Ttt'^^'^^W) = X^ and Tt{W) = 7(m,Tr(X{i)), 
where m = t{v). The first conditions shows that the above defines a splitting. To see that 
it is the one that corresponds with the vertical lift connection, we give the actions of 7^ 
and uj^ on the basis vector fields. We have Tt{X^) =_XiOT, Tt{X^) = and Tt{E^) = 0. 
Likewise, TT(Xf) = ^ o r , Tt(xY) = and Tt(E^) = 0, where r : TM/G -^ M/G. It 
follows that 

7^(Xf)=Xf, 7^(X7)=X7, 7^(^)=4^ 

and from these formulas we get for the associated left splitting 

o;^(Xf) = 0, ^^(Xn=0, u;-iE^J=E- and u^E:) = = u^E:) . 



13 



From the above relations it is clear that w^ = t*w, and that therefore 7^ is indeed the 
right splitting corresponding to the vertical lift of the principal connection on M specified 
hj zu, as it was defined initially. 

One can find a connection whose definition is somewhat similar to that of the vertical 
lift connection in [8], albeit in a much less general context: the authors deal only with a 
Lagrangian system with symmetry and restrict the dynamics to a particular value of the 
momentum map. 

The second connection of interest is a connection on the middle vertical line in the M 
square diagram, that is, it is a connection on the bundle Tvr*^ : TM — > T{M/G). It is in 
fact a particular case of a quite general construction which can be described as follows. 

We first make an obvious remark. The complete lift operation X{M) —>■ X{TM), X i— > 
X^, is not C'^{M)-\meai: in fact for a function / on M we have {fXf = /X^ + fX^, 
where / is the total derivative of /; the point to note is that {fX)'^ is a C°°(TM)-linear 
combination of X'^ and X^ . Suppose now that M is equipped with a distribution (vector 
field system) V. Let {X.i\ be a local vector field basis for X>, and consider the local vector 
fields {Xf,XJ} on TM: they are linearly independent, and there are 2dimP of them. 
Furthermore, if {Yi\ is another local basis for V then the span of {Y[^,YJ} coincides 
with the span of {Xf,XJ}, as follows from the observation above about (/X)^. The 
span of {Xf,XJ}, where {Xj} is any local basis of V, accordingly defines a 2dimX>- 
dimensional distribution V on TM . Suppose next that is a diffeomorphism of M 
and 0*^ is the induced diffeomorphism of TM. Denote by (f)'^ the action of (f) on vector 
fields on M, 0^^ the action of 0° on vector fields on TM. Then 0-^(X)^ = 0^'^(X^) 
and (j)'^{XY = 0° (^^) (these are the integrated versions of two formulas for brackets 
between complete and vertical lifts which we used earlier). Thus if V is invariant under 
the action of some group G on M then T>' is invariant under the induced action of G on 
TM. Now let M -^ M/G be a principal G-bundle and V the horizontal distribution of 
a principal connection: then V is a G-invariant distribution on TM which is transverse 
to the fibres of TM -^ T{M/G), that is, a connection on TM — > T{M/G), which is 
G-invariant in the appropriate sense. 

It is easy to describe the left splitting of the new connection, as follows. 

The complete lift construction can be extended from vector fields to tensor fields, as is 
shown in [IT]. In particular, given a type (1,1) tensor field A on a manifold M, its 
complete lift A'^ is a type (1, 1) tensor field on TM with the following properties: 

A^(X^) = A{Xy, A^(X^) = A(X)^ CxcA"" = (CxA)'', 

for any vector field X on M. Moreover, for any two type (1, 1) tensor fields A, B on M, 
A'^B'^ = {AB)^. The complete lift A'^ may be described explicitly as follows. Regard A 
as a fibre-linear map TM — > TM, fibred over the identity. Let a : TTM —>■ TTM denote 
the canonical involution: then A'^, regarded as a fibre-linear map TTM — > TTM, is given 
by A'^ = 0" o TA o a (where TA is the tangent map, or differential, of the map A). 

Proposition 3. Consider the left splitting u : TM ^ M x q of a principal connection 
on Ti^' : M — > M/G. The complete lift uj'^ of u defines a G-invariant connection on 

TM ^T{M/G). 
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Proof. As we pointed out earlier, u can be considered as a type (1, 1) tensor field on M, 
when we regard M x g as a subbundle of TM; from this point of view, for each m G M, Um 
is the projection onto the vertical subspace of TmM along the horizontal subspace; since 
it is a projection operator u satisfies u"^ = uo. The fact that the connection is principal 
is equivalent to the fact that, as a type (1, 1) tensor, oj is G-invariant, which is to say 
that uj o %l)^ = ip^ o u for all g & G, where ip'^ is the G-action on vector fields. When 
G is connected the latter condition is equivalent to C^j^jU = for all ^ G g. We take the 
complete lift u'^, to obtain a type (1, 1) tensor field on TM. Now (cu^)^ = (^^)'~^ = ^'^, 
so u'^ is a projection operator. From the formulas for the action of u'^ on vertical and 
complete lifts it is clear that it vanishes on vertical and complete lifts of vector fields 
which are horizontal with respect to u, that is, on V. Moreover, for any ^ G g we have 

^""(Cm) = {^{^mW = Cm, ^"(e.^/) = i^i^M)r = Cm, 
so that im(u;*^) can be identified with g x g in the required manner. 
Finally, we have 

At„^" = %^" = (^Cm^)' = 0, 
which expresses the G-invariance oi u^. D 



The connection determined by u'^ was first described, in its essentials, by Vilms [16]. In 
fact it was shown in [15] that a connection on a vector bundle E —>■ X induces a connection 
on the bundle TE —>■ TX. Of course the bundle vr*^ : M -^ M/G we are dealing with in 
the current situation is not a vector bundle; nevertheless, Vilms's result may be extended 
to cover it. We therefore call this connection the Vilms connection; however, we denote 
its splittings by 7° and u^ (as before). 

Note the important but somewhat subtle difference between the constructions of the two 
connections: in constructing the vertical lift connection we specify the initial connection 
by the g-valued connection form w, but in constructing the Vilms connection we specify 
it by the type (1, 1) tensor field u defining the right splitting. We mention this because 
there is a concept of the vertical lift of a type (1, 1) tensor field, and it is important to 
reahse that we do not use this concept here. 

The right splitting 7*^ of the Vilms connection is a map (T7t'^')*TT{M/G) -^ TTM, which 
may be specified as follows. We denote by a the canonical involution of TT{M/G). Let 
TTt"{v) = V. The right splitting of the Vilms connection is given by 

7°:(t;,F^)^a(T7(t;,a(F^))). 

The fact that this is a splitting is due to the property TTti"' oa = aoTTn'^', as it is easy 
to see. Indeed, 

TT7r"o7°(t;,F) = aoT{Tn" o^){v,a{Y)) = aoa{Y) =Y. 

We calculate the corresponding right splitting, and confirm that it is u'^. The right 
splitting is given by id — 7'^' o TTvr*^; we have 

id-^^oTTn"' = id-aoT-foaoTTn^' 

= id — o" o T7 o TT-K^' o o" = id — o" o T(7 o T-k^') o a 
= id — cr o T(id — uj) o a = a o Tuj o a = uj*^ 
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In terms of the standard basis {Xf , E^, X^ , E^} we have 

c.^(Xf) = 0, c.^(X7) = 0, u;^{E-) = E- and u;-{E:] 



El 



equally, cuC(E^) = E^. Since TT-n^'i^Xf 
also follows that 



9x' 



7^ 



A)=xf-^-(Xf) = Xf and im-]=Xj-u-{X^) = Xl 



dx' 



Tvr*^ and TTtt^'{X^ 



d 



i^ ° ^^"> it 



dv^ 



We noted above that because of the invariance of a;, the right splitting uj"^ of the Vilms 
connection is invariant. It follows that the Vilms connection can be quotiented to give a 
connection on the M/G-square. 

As a consequence of the existence of the two connections described so far, we can deduce 
for both square diagrams a third connection. Clearly im(7°) C im(7^), so there is a 
connection 7' such that 7° = 7^ o 7', and a connection 7' such that and 7^ = 7V o 7'. For 
the appropriate bases 



'''^^-''' 



"^d '''(^) = -V, 



and 



a;'(Xf) = 0, cu'(X7) = and ^.'(^1 



^:: 



4 Second-order systems 



We now come to the consideration of second-order systems. We assume given a second- 
order differential equation field, that is, a vector field F on TM such that TtT{v) = v 
for all V G TM, where r : TM ^ M is the tangent bundle projection. Furthermore, 
we assume that F is G-invariant, so that it satisfies ■?/'J^^'^F(f ) = F(-?/'™-'f). There is 
therefore a section F of TTM/G such that tt^^*'^ o F = F o tt^^'^. Under the appropriate 
maps F projects onto Fi and F2 as shown below, and can be decomposed into elements 
which are boxed in the diagram. Analogously, F projects on Fi and F2, and has a similar 
decomposition. Of course, all the elements are related to each other in an appropriate 
way; for example Fi G Sec{{n™yT_{TM/G)), which is given by Ti{v) = {v,Ttt^^'{T{v))), 
can be reduced to the vector field Fi G X{TM/G). 

Sections of bundles over M : Sections of bundles over M/G: 

. c.-(r) ' ' 



u;^r) 



cu'iT,) 



co-iT) 



c.-(r) 



rp„TM 

-^^- r, 



c.v(r) 



u-{T) 







[T. 



TTtt" 



{tt'^^)*Tq 



[TTit'-\ 



u'{T,) 



Ti 



Tq 
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In fact the connections give a G-invariant decomposition of F into three parts: 

r = 7"(ri)+a;^(r) 
= 7"(r2) + ^"(r) 

Since F is a second-order differential equation field the 0- valued function c^^(F) on TTM 
is given by 

tu''{T{v)) = zu{Tt{T{v)) = U7{v) 

for all V G TM; thus (for a given choice of w) tz7^(F) is the same for all second-order 
differential equation fields F. 

There is an analogous three-way decomposition of F, as shown in the right-hand diagram. 
Here Fi is a vector field on TM/G, F2 is a section of q*TT{M/G) (i.e. a vector field along 
q) and ^(Fi) G Sec (TM/G x g). 

As well as being a section of TTM/G -^ TM/G, F may also be regarded as a section of 
the so-called prolongation bundle T^{TM/G) -^ TM/G, whose fibre at -u G TM/G is 

TliTM/G) = {{w,Xi^) G TM/G x T{TM/G) \ q{w) = Tt{X^)}. 

Theorem 9.1 of j6] shows that the identification of the quotient bundle TTM/G —>■ TM/G 
with the above bundle is given by the isomorphism 

TTM/G -^ T^{TM/G), [Wj ^ {7f™TT{W),T7c™{W)). 

This map is independent of the choice of W since for any other ip^'^'^'W within the same 
equivalence class, Tr^ipg^'^'W) = T{t <^ ^J™ ){W) = Tt{W). The second component of 
the isomorphism is in fact the map [Tvr^*^] : TTM/G -^ T{TM/G). Keeping in mind 
that F is a second-order differential equation field, F G Sec(TTM/G) can be identified 
with the section v h- *■ {v,Ti{v)) of the prolongation bundle. As a consequence, the com- 
posing parts of this section satisfy Tt{Vi{v)) = g{v). This property clearly resembles 
the defining property Tr(F(f)) = f of a second-order differential equation field; sections 
of the prolongation bundle of the above form were therefore called 'pseudo second-order 
differential equation sections' in e.g. [12] or 'second-order differential equations' in e.g. [6]. 

Next, we will discuss the reconstruction process. We will use the following notations. 
Let v(t) G TM denote an integral curve of F and let c(t) be the corresponding base 
integral curve, that is, c{t) = r(f (t)) G M. It follows from the fact that F is a second- 
order differential equation field that v = c (when we consider the latter as a curve in 
TM). We will write v{t) = 7r^*^(t;(t)) G TM/G and c{t) = 7r*^(c(t)) G M/G. Obviously, 
t{v) = c and moreover T'K^"{v{t)) = g{v{t)) = c{t) G T{M/G). In a previous section 
we encountered the horizontal lift of c (with respect to the connection on vr*^), which we 
denoted by cP. 

We first note that the vertical lift connection is a principal connection on the principal 
fibre bundle tt^*^ : TM — >■ TM/G. So just as in the first-order case we can construct the 
integral curve t 1-^ v{t) of the invariant vector field F G X{TM) from an integral curve 
t ^ v{t) = TT™{v{t)) of the reduced vector field Fi G X{TM/G), by 
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• taking the horizontal hft v^ of v through f (0) (with respect to 7^), and 

• finding the solution t ^^ g{t) E G oi the equation 

e{g) = w''{Tovr) (11) 

with (yf(0) = e (where 9 is the Maurer-Cartan form of G); 
the required integral curve is given by 

The right-hand side of equation ( ITTl) can equally well be written as w{v'^). 

Let us look at the relation that determines the horizontal lift v"' of v to TM. Let w be a 
given curve in TM/G (not necessarily an integral curve of Fi). By definition, v"' projects 
onto V and is a solution of 

From the conditions that determine 7^, we know that this is equivalent with the properties 
Ttc'^"ov = V, and Troy = 7(Tro{;). The first property simply recalls that 7r^"({;'^) = v. 
If we denote as before t o v = c, then Tt o v = c. So, we can deduce from the second 
property that Tt o v = c . To conclude, the curve v'' is completely determined by the 
properties t o v'' = c^ and vr^" o v'^ = v. 

Any element of the vector bundle TM/G -^ MjG can be written as a sum of two parts 
via the splitting 7 of the Atiyah sequence ([2]): 

^ = ? + 7(c), 

where ^ = '0{y) G g. We can use this to give a more explicit formulation of iP . Let 
^(t) G be such that ^ = [cf,.^]. Then the curve ^m o c"^ + c in TM projects onto cP" 
by means of r and projects onto v = ^ ^ 7(c) by means of tt"^*^; so it can only be v'^ . 
Therefore, w{y'~^^ = ^ and the horizontal part of 'u''' is c . 

Suppose now again that v{t) is an integral curve of Fi: 

v = Tiod; (12) 

and that v(t) = ''P'^u^v^it) is an integral curve of F. Notice that 

that is to say, the curve in G required to bring v"^ to v in TM is the same as the curve in G 
required to bring cP" to c = r oti in M. In fact from equation (jl]), since F is a second-order 
differential equation field, 

v = c = r/'mg))MO-c' + f). (13) 

So, in the case that v is an integral curve of Fi, the curve ^{t) = w{v'^(t)) G g must equal 
0{g), which agrees with equation (ITTi) . 
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We turn finally to the integral curves v oiTi. We will use the connection u' to decompose 
equation flT2l) into two coupled equations for the two curves ^ G and c G M/G that 
constitute v. The first equation is related to r2 G Sec{g*TT{M/G)), which can be 
considered as a map T2 : TM/G -^ TT{M/G); thus fa o {) is a curve in TT{M/G). 
The second equation is related to a;'(ri), which is a vertical vector field on TM/G. This 
vector field can be regarded as a section of (r)*g, thus as a map TM/G — ^ g; so ci;'(ri ov) 
can be regarded as a curve in g. The projection of this curve onto M/G is obviously c. 

We need to introduce one more concept: that of the associated linear connection on the 
associated bundle q -^ M/G (see also [H [7]). In fact, we will only need its covariant 
derivative operator -^ which acts on curves ^ in q. Let c(t) be the projection of ^(t) on 
M/G and let c(t) be any curve in M that projects on c{t). Let ^(t) G g be such that 
C, = |c, ^]. Then, the covariant derivative of ^ can be defined as 



Dt 



|c, ^ - [w o c, . 



(^ stands here for the projection on the second argument of this curve in Tq = x g). To 
see that this definition is independent of the choice of the representative in the equivalence 
class, take any other d{t) G M with d{t) = '4^h(t)'^i^) ■ "^^^ corresponding curve ^'^(t) in 
such that ^ = {d,^'^] is then equal to adh-i^. Moreover, d = -ip™ [c + {9{h)){c)) and 
C,d = adh'i (^ + [6{h),Q). So, indeed, 

[d, ^^ -[mod, ell = Id, adu-. [i + [e{h),^]) - [a4-i (w o c + e{h)) , a^-i^]! 
= K^c, a4-i (^ - [^ o c, e])] = [c,i-[rxjo c, e]]. 

Remark that in the particular case of the horizontal lift, '^[c'^,^] = [c''^)^]- 
One can show that the explicit formula for the associated linear connection is 

V^ : X{M/G) X Sec(0) ^ Sec(0) : {X,l) ^ V^l = [l(X),^]. 

Here [■, ■] stands for the above mentioned Lie algebroid bracket of the Atiyah algebroid 
TM/G. For more details, see e.g. [7]. 

Theorem 4. Let ^{t) G 0, c(t) G M/G and put v = C, + 7(c). Ifc and ^ are solutions of 

C = T2OV, 

D^l -,- , (14) 

-m = "'^^^°")' 

then V is an integral curve o/Fi. Solve c = 7(0^,0) for c^{t) G M and let ^(t) E Q he 
such that e = |c^,e]- If g{t) & G is a solution of 

0{9) = e (15) 

then the curve v = ip™{S,M o c'^ + c ) = ip^v"' is an integral curve ofV. 

Conversely, suppose that v is an integral curve ofV. Let v = vr^" o v, c = n'^' or ov and 
^ = uJ o V. Then v is an integral curve ofTi and c and ^ satisfy [T4\ )- Compute cf from 
c = 7(0^,0). Let g E G be such that c = ip^'c' and let ^ E g be such that ^ = lc',Q. 
Then g satisfies equation ( 173]) . 
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If (i is a curve in M such that TT^'{d) = c and if we denote by C,'^ the curve in q which 
is such that ^ = {d, ^] , then the last equation of (fT^ could equivalently be written as 
^ — [tz7((i),^] = u'(Ti o v'^), where v'^ = C,m ° d + 'y{d,c) is the unique curve on TM that 
projects on both v and d. Indeed, the relation between uj'{Ti) and uj'{Ti) is 

and therefore, cJ'(ri o v) = {d, u'iTi o v'^)]. 

If one is interested only in the coordinates on M/G ('shape' variables in [2]), it is necessary 



only to solve equations flT4j) where the symmetry has already been cancelled out. If the 
whole motion on M is required one will have to solve the whole system. 

The proof of theorem will follow from the considerations of the coordinate version of the 
reduced second-order equations in the following paragraphs. 

We begin our description of the coordinate expression of the equations with a general 
remark. If we take any local basis {Xa} of vector fields on some manifold M, not nec- 
essarily a coordinate basis, and express any tangent vector v at m G M in terms of this 
basis so that v = v°'Xa\m, then the f" will serve as fibre coordinates on TM. A vector 
field r on TM will be a second-order differential equation field if and only if it takes the 
form r = v"X^ + F^X^, and its integral curves will satisfy v" = F". In terms of a 
new basis {Ya}, where Ya = A^^Xp, we have v = v^Xalm = w°'Ya\m where w'^A'^ = v". 
Moreover, F = w°'Y^ + C^Y^ where F" = A'aG^ + A'^w'^, the overdot here indicating the 
total derivative. 

We turn now to the case of interest. We have defined on M two local vector field bases 
{Xi, Ea} and {Xi, Ea}, with 

where the x* are local coordinates on M/G, and {Ea} is a basis for q. Both {Ea} and 
{Ea} are bases of vector fields which are vertical with respect to the projection ii'^' : 
M -^ M/G. The first, which consists of fundamental vector fields, we called the moving 
basis. The body-fixed local basis {Ea}, on the other hand, consists of G-invariant vertical 
vector fields. We have Ea = A\Ei, where the coefficients A\ satisfy Ea{Al) + C^rfA^ = 
(equation ([6])). 

For any v G TmM we set 

V = v'Xi\m + V^'EaU = v'Xi\^ + W^'EaW, V' = A>^ 

The w* may be regarded as the fibre coordinates on T{M/G) corresponding to the base 
coordinates x*. We show in the following paragraph that the f" satisfy 

E^{v'^) + CSy = 0; (16) 

so from equation ( TTOj) we may consider v i-^ v"'Ea as defining a section of t*q -^ TM/G; 
this is just the section u^(T). On the other hand, since v"" = A'^w^, where the A'^ are 
functions on M , we have 

\a T?Ct „,,c 



AlE^,{w^)^E,{Al)w^^Cls^^ = ^- 
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but from equation (Q 

whence E^{w°') = 0, as one might have expected. 

Equation fITBl) is a consequence of the following general considerations. Let {Za} be any 
local basis of vector fields on a manifold M, with dual basis of 1-forms 6°. Let 6 be 
the fibre-linear function on TM defined by a 1-form 6 on M, so that 6{x,v) = 9x{v). 
Then if t;" = ^°(a;,t;), v = v"Za\x- For any vector field Z on M, Z^(^) = C^O. But 
Cz^ef^iZy) = -e^i[Z^, Z^]), so Cz-0^ = -C^^O^ where [Z^, Z^] = C^^Zp. That is to say, 
Z^{v^) = -Cl^^v^. It follows that [Z^,v'^Z^] = 0. 

The second-order differential equation field F may be written 

r = v'X^ + v^E^ + D'Xj + D''E^ 
= v'Xf + w^'E': + D'X7 + F"E 



a ) 



we have 



^a ^ j^apb ^ j^a b_ 



By assumption F is G-invariant, which is to say that [-^ajF] = 0. Now [E^^v'^Xf] = 0, 
and [E^,v^E^] = 0, as follows from equation (ITB]) and the argument that establishes 
it. Moreover [E^,w''E^] = since E^{w'') = and [E^,E^] = 0. Thus v^E^ and 
w'^E^ are both G-invariant; they are not however equal, but differ by the vertical vector 
field w'^A'^E'^, which accordingly is G-invariant. Next, [E^,D^X^] = E^{D'^)X^ since 
[E^,X^] = [^a,Xi]^ = 0. On the other hand, 

[E-,D'E^] = E-{D')E^ + D'[E-,E^] 

= [mD') + D^Cl) El 
Finally, [E^,E^El = E^{E^)El The remaining coefficients of F must therefore satisfy 

that is to say, D^ and F^ are G-invariant, while the D"^ may be regarded as the components 
of a section of T*g -^ TM/G. 

Observe that 

cu^(F) = v^E^ e X{TM). 

The corresponding g-valued function vo^iV) is given by w^(F) = v'^Ea, it is independent 
of the choice of F, as we remarked before, and as we showed above it in fact determines 
a section of r*g. 

We may also express F in terms of the mixed basis: 

F = v'Xf + v^E"^ + D'XJ + CEl, 
where 
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the A^ being the components of the inverse of the matrix {A'^). It is easy to see that 

E^iG'') = 0. We have 

Ti = v'X^ + D'X^ + CE^. 

Then cu'(ri) = G^E^ G Sec{{n^^'yT{TM/G)), and so 

7V(cu'(ri)) = G"i; = D^'E^ e X(TM). 

The three-way decomposition of F at the level of the M square diagram is therefore given 
by 

r = 7"(r2) + 7"(^'(ri)) + a;^(r) 

= {v'Xf + D'Xl) + D^El + v^'E'^. 
Among the equations for the integral curves of F we find 

x' = D\ 



w" + A'iAiw'' = G". 



We can write the latter as 

Since the w" are G-invariant they can be taken, together with x' and f', as coordinates 
on TM/G. The integral curves of Fi are the solutions of the equations 

x' = D\ 



w'^ + AM^^ = G"^. 



This latter equation has a familiar structure: one could think of the term A'^A'^ as repre- 
senting the 'angular velocity' of the body-fixed frame with respect to the moving frame, 
and w"' as components of some velocity with respect to the body-fixed frame; the whole 
term A'^A'^w'' is then of Coriolis type. 

We can also write the same equation as 

A^w'' + A^w^ = AIG\ 

which is equivalent to i)"" = D"^. 

Finally, if we rewrite Al as x^ Xj{Al) + v'''Ed{Al) and use the formulae Xj{A'^ = •y'^C^bAl 
obtained earlier (equation ([7])) and Edi^A"^) = —C^^Al (equation ([6])) we find that 

AlAlw' = A:{x^jP^,,A:-v'C^,,AI)w' 

= x^jfC^,,w' - Ayv^C^,, = x^jfC-,,w\ 



But ■jjG^fj = Tj*^ are the connection coefficients of the adjoint connection; the equation 
for w"' is therefore equivalent to 

w'' + T^'f^x'w^ = G^ 

which is in agreement with the second of equations (1141) in the theorem. 
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5 An example 

In this final section we determine the reduced equations for an interesting class of second- 
order differential equation fields. 

The case to be discussed is that in which there is a 'kinetic energy' metric k on M, with 
Levi-Civita covariant derivative V, and the equations of motion of the original dynamical 
system take the form 

VcC = F(c, c) 

for the curve t i— *> c(t) on M. Such a system may be called a system of mechanical type, 
with F representing a force field. We hasten to point out, however, that according to the 
philosophy of the paper as we described it in the Introduction these features of the system 
are incidental to our main purpose, which is to illustrate the methods described above 
using a familiar example, rather than to discover properties of systems of mechanical type 
that they have because they are systems of mechanical type. 

There is great potential for confusion here, since we will now have two connections of 
fundamental importance to deal with, the Levi-Civita connection and the connection on 
the principal bundle vr*^ : M -^ M/G (when we have defined the group G and its action); 
we warn the reader to be on guard. 

We may write F in the form F = Fq + $ where Fq is the geodesic spray of the Levi-Civita 
connection and $ is the force term on the right-hand side of the equations of motion 
considered as a vertical vector field on TM. We now examine the possible symmetry 
conditions. For any vector field Z on M and any affine spray Fq, [Z'^,Fo] is vertical and 
quadratic in the fibre coordinates. On the other hand, [Z'^,*!)] is vertical since $ is; but 
in cases of interest (for example, when F is independent of velocities, or linear in them, 
or a combination of the two) there will be no terms quadratic in the fibre coordinates; so 
it is natural to consider the situation where [Z'^,Fo] and [Z^,<I>] vanish separately. Now 
[Z^, Fq] vanishes if and only if Z is an infinitesimal affine transformation of the symmetric 
covariant derivative defined by Fq, which in the case under discussion is the Levi-Civita 
connection of k. Since any infinitesimal isometry is affine, it is natural to assume further 
that G is a group of isometrics of k, whose elements in addition leave invariant the force 
term, as represented by the vertical vector field $. Such a group is always a symmetry 
group of F, and in many cases the maximal symmetry group will be of this form. 

We now turn to the choice of a vector field basis on M adapted to the group action. In 
this case there is a natural choice for the connection on tt'^': take its horizontal subspaces 
to be the orthogonal complements of the tangent planes to the group orbits; they are 
G-invariant since the group consists of isometrics. The vertical vector fields Ea comprise 
a basis for the Killing fields or infinitesimal isometrics. We shall however work with an 
invariant, body-fixed basis for the vertical vector fields; that is, we choose a local basis 
of vector fields of the form {Xj, Ea}. The components of k in this basis are denoted by 
kab, kai, kij in the obvious fashion. The kab are evidently G-invariant. By construction, 
kai = 0. The kij are also G- invariant, and so define functions kij on M/G which are 
the components of the reduced metric, say k, with respect to the local vector field basis 
there. We may without loss of generality take this basis to consist of coordinate fields, as 
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before; the Xj will not in general commute, but [Xj,Xj] will have components tangent to 
the group orbits; we set [Xj, Xj] = K^^Ea (this in effect defines K as the curvature of the 
connection) . 

The connection on tt" has now been entirely taken care of; references to a connection 
henceforth always mean the Levi-Civita connection. 

We set $ = $*X7 + "^"-EJ; by assumption both $* and $" are G-invariant. 

In order to find the reduced system it is necessary to express V in terms of the adapted 
basis. For this purpose we need the Christoffel symbols of the Levi-Civita connection 
with respect to the basis {Xj,^^}: we set 

and so on. The order of indices is important; though the Levi-Civita connection is sym- 
metric, it is represented here with respect to a non-coordinate frame. To calculate the 
Christoffel symbols we need the brackets of the basis vector fields. Recall from equation ([8]) 
that [Ea.Xi] = •y^C^i^Ec, and that the connection coefficients of the adjoint connection 
are given by T^^ = jfC^!^. We therefore have the following bracket relations: 

[Ea, Eb] = -C^f^Ec, [Xi, Ea] = TiaEb, [Xi, Xj] = K^jEa. 

Since all of the vector fields appearing are G-invariant, so are all of the coefficients on 
the right-hand sides. Furthermore, since all of the brackets are vertical the Christoffel 
symbols with upper index i will be symmetric in their lower indices. 

Using these data in the standard Koszul formulae for the Levi-Civita connection coeffi- 
cients of k with respect to the basis {Ea,^i} we find that 

be ~ 2 (~^bc + ^ \f^be(^dc~^ ^ceC'bd)} 

r^, = ik'^{-x,{hc) + kMrf, + h,Tf,) 

Fj^ = iA;-(X,(A;b,)-A;,,T;i + fc,,T;i) 

F,^. = lk'''{x,{kbc)-kurf,-k,drf^ 

jb ~ 2^ '^bc-l^jk — '-bj 

pa _ i /<"" 
jk 2 jk 

-pi _ p« 
jk '-jk^ 

where in the final line the T,^ are the Christoffel symbols of the reduced metric kij. 
It follows that 

Fo = x^Xf + w'^i^ 

- {x^^V;, + x^n^iV;, + Vi^)+ w'w^Vi) XJ 

- {x^x'T^, + x^w\T^, + F,^.) + w'w'^T^:) Ej 
= x^X^ + w^'E^ 

yx%i - x^w'k^'kbcK]^ + w''w'k'^ {-IXjihe) + hdTf^)^ XJ 

- [x^w'k'^' [X,{h,) - kuTf,) + w^w'k'^^h.Ci;) El. 
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The reduced equations are therefore 

x' + Vlfj:^i^ = $* + x^w''k"'hcK^k + w'w'k'' ilXjihc) - hdTf^) 
w'' + T;1,x^w' = •I.'^ - x^w'k'^'^ (X,(fc,J - k,,Tf^ - fcedTji) - w'w'k''%,Cl. 

The first equation can be written 

where D/Dt is the covariant derivative operator of the Levi-Civita connection of k, and 
$j = kij^^ . We can write the equation for w" in either of the following two forms: 

Dt Dt '^^ 

^ = <^^-w,w^k^'^Cl, 

to obtain the second we have used kab to lower indices. 

When F = 0, that is, when $* = $" = 0, we obtain Wong's equations jUIll]. The case in 
which F 7^ but G is l-dimensional is discussed by Bullo and Lewis ^. We shall show 
that in both cases our equations subsume those of the cited authors. 

In the case discussed in [H [H] , in addition to F = it is assumed that the vertical part of 
the metric comes from a bi-invariant metric on the Lie group G. This means in the first 
place that £^ k{Eai Eh) = as well as C^ k{Ea, E^) = 0, and secondly that the kab must 
be independent of the x*. From the first condition we easily find that the kab must satisfy 
kadG^^ + kbdCac = 0, and therefore fcacTTff, + kbc'^fa = 0- From both together we see that 
the kab must be constants. Thus ^g/g^iikbc) = and w^w^kbeG^^ = —w'^w^kdeGl^ = 0, 
and the reduced equations are 

D(kijx^) 



Dt 

D^Wa 



-x'wbKt^ 



0. 



Dt 
These are equivalent to the equations given in [U [H]. 

In the l-dimensional case we have a single Killing field F; this vector field is also clearly 
invariant, so we shall simplify the notation by denoting it simply by F (Bullo and Lewis in 
fact write X for this vector field). There is but one component of kab-, which is fc(F, F) = 
|Fp, and k with notional upper indices is just |F|~^. Furthermore, |Fp is itself invariant, 
and may therefore be considered as a function on M/G. An arbitrary tangent vector 
V may be written in the form V = vE + f *Xj (so v is to be identified with the single 
component of w"-), and since the Xj are orthogonal to F we have 

k{V,E) 

'=^Er'^ 
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Bullo and Lewis call the map V ^— *> k{V,E) the momentum map and denote it by Jx- 
They also introduce a type (1, 1) tensor field on M/G which they call the gyroscopic 
tensor, which they denote by Cx- The gyroscopic tensor is given essentially as follows. 
The covariant differential VE is a type (1, 1) tensor field on M. Let us denote by E-'- 
the distribution orthogonal to E, that is, the distribution spanned by the vector fields 
Xj. Then VE may be used to define an operator on E-^, by first restricting its arguments 
to lie in this distribution, and then perpendicularly projecting its values into it. Now in 
general we have 

Vx.Ea = T^^Xj + Ti^Eb] 

so we are concerned here with F/^ = —^kabk^'^K'^i^, albeit in the 1-dimensional case. In 
fact if we write [Xj,Xj] = KijE the gyroscopic tensor in component form is 

ci = lEi^k^'K.k = i^rFV^fc. 

It is clear from this that Cj is invariant and that Cij = k-ikC^ is skew- symmetric. Moreover 
in the 1-dimensional case T^^ = 0. The reduced equations of motion in this case are 
therefore 

x' + t'x^x^ = ¥ + vCix^ + WV'^]^ 

11 = <J> r — ■ — — 

\E\^ dx^ ■ 

Here ^^ is the ii^-component of the force. These equations agree with those given by Bullo 
and Lewis. These authors deal mainly with the case in which the force is derived from a 
potential, and the last term on the right-hand side of the first equation is subsumed by 
them into the so-called effective potential. Bullo and Lewis actually give two versions of 
the reduced equations: one is in terms of v, and is derived above; the other is in terms 
of /i = l-Epf , and the second of the reduced equations is then simply /i = 0. But since 
lE'P is the single component of kab, /i is just the single component of Wa- Furthermore 
Bullo and Lewis have $„ = 0. So the equation /i = is just the second reduced equation 
written in terms of Wa- 

There is a simple explicit example which nicely illustrates both of these cases, namely the 
Kaluza-Klein formulation of the equations of motion of a charged particle in a magnetic 
field. The Hamiltonian and Lagrangian approaches to this topic are well-known: see for 
example [9]; here we derive the equations from those given above. For M we take E^ x S, 
with coordinates {x'^,6). Let Ai be the components of a covector field on E^, and define 
a metric k on M, the Kaluza-Klein metric, by 

k = 5ijdx' dx^ + {Aidx' + dOf 

where {5ij) is the Euclidean metric. The Kaluza-Klein metric admits the Killing field 
E = d/dd. The vector fields Xi = d/dx^ — Aid/dd are orthogonal to E and invariant; 
moreover kij = k{Xi,Xj) = 6ij, while \E\ = 1. Finally 



[Xi,Xj] 



dAi dAj\ d 



dx^ dx^ I do 
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Putting these values into the reduced equations above we obtain 

X = VX-' — — — — ^ , V = \}. 

These are the equations of motion of a particle of unit mass and charge w in a magnetic 
field whose vector potential is Aidx^. 
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